Anderson impurity model in a semiconductor 
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We consider an Anderson impurity model in which the locally correlated orbital is coupled to a 
host with a gapped density of states. Single-particle dynamics are studied, within a perturbative 
framework that includes both explicit second-order perturbation theory and self-consistent pertur- 
bation theory to all orders in the interaction. Away from particle-hole symmetry the system is 
shown to be a generalized Fermi liquid (GFL) in the sense of being perturbatively connectable to 
the non-interacting limit; and the exact Friedel sum rule for the GFL phase is obtained. We show 
by contrast that the particle-hole symmetric point of the model is not perturbatively connected to 
the non- interacting limit, and as such is a non- Fermi liquid for all non-zero gaps. Our conclusions 
are in agreement with NRG studies of the problem. 



I. INTRODUCTION 

The Anderson impurity model (AIM), in which a 
single, locally correlated orbital couples to a non- 
interacting metallic band of electrons, is a longstanding 
paradigm of strongly-correlated electron physics. Con- 
ceived originally^ to explain the formation of localized 
magnetic moments on impurities in non-magnetic hosts, 
it has since formed the cornerstone of our understanding 
of the Kondo effect* and related many-body phenomena. 
Interest in the area is currently particularly strong, both 
experimentally and theoretically, after the Kondo effect 
was predicted 3 -! 4 - and then directly confirmed 5 -* 6 - to arise 
in mesoscopic quantum dot systems.— 

After some 50 years of intense theoretical work, the 
spin-i Kondo effect as manifest in Anderson's original 
model is naturally rather well understood?. Below some 
characteristic Kondo temperature Tk, a complex many- 
body state develops in which the impurity spin is com- 
pletely screened by the host metal, leading at low energies 
to a 'local' Fermi-liquid and universal transport proper- 
ties. 

Being a low-energy phenomenon, the Kondo effect is 
of course crucially dependent on both conduction band 
states near the Fermi level and the low-energy spin de- 
grees of freedom of the impurity. This has inspired much 
research into other quantum impurity models — involving 
more complex impurities and/or host densities of states — 
with the aim of identifying the various types of Kondo 
effect that may arise, the conditions under which they 
do so, and the novel physics that results when Kondo 
screening cannot be achieved^. 

Here we consider the notionally simple problem of an 
Anderson impurity in a gapped host, where the density 
of states vanishes over a finite range about the chemi- 
cal potential, a model not only of relevance to Anderson 
impurities in semiconductors but alsoi - to the topical 
issue of impurities in BCS superconductor s n i 12 . In re- 
moving the all- important low- lying states of the host, one 
would certainly expect the Kondo effect to be precluded 
for large enough gaps: the question is, can the effect still 
arise for sufficiently-small gaps, or is it destroyed as soon 
as a gap is opened? 



This question has indeed been the subject of a num- 
ber of previous papers. Poor man's scaling, the 1/iV ex- 
pansion and the non-crossing approximation predic t 13 ! 14 
that the Kondo effect always arises whenever the gap is 
less than the Kondo temperature in the absence of the 
gap, while for larger gaps the system undergoes a quan- 
tum phase transition to an 'local moment' (LM) phase 
where the impurity spin remains unscreened as T — > 0. 
In addition the problem has been studied numerically 
by the density-matrix renormalization group 15 and quan- 
tum Monte Carlo 13 i 14 ' 16 , but with no general consensus 
reached regarding the nature of the quantum phase tran- 
sition. The numerical renormalization group (NRG) 16 ' 17 
on the other hand has been used to argue that the Fermi- 
liquid regime associated with the Kondo effect exists only 
away from particle hole-symmetry, and then only below 
a certain critical gap. In the particle-hole-symmetric 
limit it is foun d 16 ' 17 that the Kondo effect never arises 
and the ground state is the doubly-degenerate LM phase 
for arbitrarily small gaps. 

In this paper we study the problem analytically, within 
a perturbative framework which includes both explicit 
second-order perturbation theory and self-consistent per- 
turbation theory to all orders a la Luttinger— . In ad- 
dition to confirming the basic predictions of the NRG 
study-i 7 -, our analysis provides a number of exact re- 
sults, including the analogue of the Friedel sum rule, 
which serve as touchstones for approximate theories of 
the gapped AIM (GAIM). In a subsequent paper—, we 
present a local moment approach^ - to the problem, the 
results of which agree very well with the conclusions of 
the present work. 



II. PRELIMINARIES 

In standard notation the generic Anderson 
Hamiltonian 2 is 

H = Xl ek " kff+ 5Z eifliCT+c/ " i T"n+Xl^k c L c i ( T+ h - c -) 



where hj C 



= A 



(l) 

c- ff is the number operator for cr-spin 



electrons on 'site' j (with j = i referring to the im- 



2 



purity site and j = k to the host band states). The 
first term in eqn. |T]) thus describes the non-interacting 
host band, the second and third terms describe the impu- 
rity with onsite Coulomb interaction U , and the fourth 
term hybridises the two. For a symmetric host band, the 
particle-hole symmetric limit corresponds to the special 
point £i = —U/2 (where H is invariant under a particle- 
hole transformation). 

The dynamics of the model will be obtained from the 
retarded Green function Gij. a (uj) 



G ij . ! „(t) = -ie(t)({c ia (t),cl}), 



(2) 



differentiation of which leads straightforwardly to its 
equations-of-motion 2 ; from which the impurity-diagonal 
Green function in the non- interacting (U = 0) limit fol- 
lows. Its Fourier transform, denoted by ga- <<7 (u)), is 



1 



with A(ui) the host-impurity hybridisation function 



A(«) = 



E 

k 



ik 



£k 



(3) 



(4) 



and z = uj + b] with r\ = + a positive infinitesimal. The 
'full' and non-interacting Green functions are related in 
the usual way by Dyson's equation 



Gii ;CT (w) 1 = ffii;<r(w) 1 - £ ii;cr (w) 



(5) 



with Eii ;tT (w) the conventional (single) self-energy. 

It is convenient below to exploit the analytic structure 
of the impurity Green functions and their constituent 
parts. Let F(z) be analytic on the real axis and in the 
upper half plane, tending to zero as \z\ — > 00. Then, 
with F(ui) = F r (lu) — iF l (u>) for real lu, one has the well 
known dispersion relation 



F R (cu) 



(6) 



(with P denoting a principal value), and the spectral rep- 
resentation 



1 r 

F(fj) = - \ 



— lu' + i?/ 



(7) 



In particular, the full impurity Green function can be 
determined entirely from its spectral function, 



MkctM = -- IniGji a (uj). 
7T 



(8) 



The results above are valid for whatever form the host 
takes in eqn. ([T]): the details of the host bandstructure 
affect only the hybridisation function A(ui). Assuming 
for simplicity that VJk = V is fixed, and writing A(lu) = 
A R (w) - iAi(w), eqn. (g]) gives 



Ai(u) = nV 2 - £k) = ttV 2 p(oj) 



(9) 
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FIG. 1: Real part of the hybridization function, Ar(w), vs. 
ui/5 (solid line). The dashed line illustrates the solution of 
eqn. (|13[) as described in the text. 



with p(lu) the host density of states. Hence, taking p(lo) 
to be constant except for a gap of full width 25 centred 
around lu — 0, we write Ai(oj) for the gapped Anderson 
model as 



Ai(w) = 



for 



< 5 



for \oj\ > S. 



(10) 



The corresponding Ar(w) (= -Ar(-w)) follows from 
the Hilbert transform eqn. © as 



Ar(w) 



A 



In 



-5 



(11) 
±6 arising 



The logarithmic divergences of Ar(w) at u) 
from the gap are shown in fig. [T] 

The explicit form of the non-interacting spectrum 
du;<r{w) is readily obtained from eqns. ((3]), (fTO]) and 
pip . Two distinct contributions arise: a continuum (or 
'band') part (^.^(uj) arising when At(oj) is non-zero 



ttAqcIIJuj) 



5)A 2 



j+S 



1-5 



(12) 



and a single pole inside the gap. This pole arises for all 
£i when 8 > 0, occuring at a frequency ujq determined by 
solution of 



- e\ = Ar(wo). 



(13) 



That this equation always has one, and only one, solution 
is guaranteed by the monotonic decrease of Ar(w) from 
00 to —00 for —5 < lu < 5, as seen from the construction 
in fig. [TJ Moreover, since Ar(0) = 0, it is clear that ujo «s 
for t\ ^ 0, respectively. In the particle-hole symmetric 
limit, where t\ = 0, the pole lies precisely at u>q = 0. 
We shall see that this zero-frequency pole is the basic 
reason for qualitatively distinct physics at particle-hole 
symmetry, compared to that elsewhere. 
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FIG. 2: Self-energy diagrams up to second order in the inter- 
action U. The solid lines represent the non-interacting prop- 
agator, and the wavy lines represent U. 



The weight of the pole, Qo, follows straightforwardly 
from a Taylor expansion of Ar(o>) as 



1 



2An 



(14) 



7T S 2 — LUq 

and hence the total non-interacting spectrum is given by 

d ii;er {uj) = d\. a (u)) + Q q 6(lv - uj ). (15) 

This non-interacting Green function of course forms 
the basis for perturbation theory in U. We begin by 
analysing the problem to second-order in U, which is 
sufficient to highlight the essential features that lead to 
non-analyticities in the particle-hole symmetric limit. 



III. SECOND-ORDER PERTURBATION 
THEORY 

To second-order in U, the self-energy S(w) (= E R (a->) — 
iE 1 ^)) is expanded diagramatically as shown in fig. [2] 
To simplify the notation, we henceforth drop the ii sub- 
scripts, since we focus on the local impurity Green func- 
tion. We also drop the spin index, as the Green functions 
obtained are independent of a. 

Translating the diagrams in fig. [2] using the Feynman 
rules for the time-ordered Green function gives 



E(w) = Wn + E 2 (w) 



(16) 



where no — 2 f_ duj d(u>) is the impurity charge in the 
non-interacting limit, 

Tj2 j-OO 

M") = 7T ^ ,0 Il(cj')g(u, + u'), (17) 

and °II(a;) is the polarisation 'bubble' 
i f°° 

°U(lj) = — / <^'g{u>')g{u>> - w). (18) 
2tt J-oo 

Using the spectral representation of the Green function 
(see eqn. (0), Im °n(c<j) is expressible entirely in terms 
of d{uj), viz. 



i /-M 

-Im °U(uj) = / du'd(oj')d(u>' 

n Jo 



M). (19) 



Likewise, E 2 (w) can be written as 
/■M 

E 2 (u;) = U 2 dw'Im °n(w') d(u - sgn(w)u/)- (20) 
Jo 

From eqns. (fT9|) and (f20|) . which hold also for the retarded 
Green function under consideration here, the second- 
order self-energy E 2 (w) can be constructed from the non- 
interacting spectrum d(uj) alone. These integrals can of 
course be evaluated numerically for any given e; and 6 to 
obtain the full frequency dependence of E 2 (w). This is 
not our aim here; instead, we focus on the low-frequency 
behaviour (which can be determined analytically), since 
it is this that contains the key physics of the problem. 



A. Away from particle-hole symmetry 

The generic behaviour away from the particle-hole 
symmetric limit can be examined by focussing on some 
given ei < 0; results for e; > follow from a straightfor- 
ward particle-hole transformation. 

Given that the pole in d(uj) lies at a frequency ujq < 
for e ; < (see sec.Ql]), it is readily shown from eqns. (fT5)) 
and (H]) that 



1 



Im °n(w) = 0(|w| + uj ) Qo d b (\uj\ + uj ) 

+ [ H 6u'd h (w')d h (Lo' - M). (21) 



Since d h {u>) has no spectral weight for \tu\ < S, the second 
term of eqn. (p?Tj) is zero for < 26. The additional 
contribution from the first term of eqn. (|2ip also contains 
a gap for \lu\ < 5 + \luo\, & n d since \u>o\ < 5 we can write 



Im °TL(uj) = 



|w| < S+ \lu \. 



(22) 



Using this in eqn. (|20[) . similar arguments can be made 
for the low- frequency form of E 2 (w); it too has a gap 
around lo — 0, given by 



E 2 H = 



-d-2\u \<LJ<2d+\u} \. (23) 



The gap in E 2 (w) enables one to deduce the salient 
behaviour of E R (w). By writing E 2 (w) = 0(\ui\ — <5)E 2 (w), 
the Hilbert transform eqn. © gives 



1 f°° 
n J -co 



,9{\iJ\-S)V\(iJ) 



(24) 



and hence 



UJ 2 - LUl 



TT 



doj' 



y\-s) 



(Ul - UJ')(UJ2 - LO 1 ) 



E 2 (c/) (25) 



If both lo\ and L02 arc within the gap, the term in square 
brackets above is non-negative for all lo' . Since E 2 (lj) is 
necessarily non-negative, it follows that 



E£(wi) > El^a) 



-6 < u>i < ui2 < S. 



(26) 
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This monotonicity is important in arguments below. 

Finally, the low-energy behaviour of X 2 (w) can be 
inserted into the Dyson equation, eqn. (|5|), to obtain 
the impurity Green function to second-order in U. For 

M<*, 



G 2 (w) 



M<<5 



1 



Wn - A R (w) - 



(27) 



and, by using similar arguments to those in the non- 
interacting limit, the monotonicity of in eqn. (|26|) 
guarantees that G 2 (w) has a single pole at a frequency 
uj p2 given by 



1^2 



±t/n = A R (w p2 ) + E^(o; p2 ). (28) 



Inside the gap the second-order Green function G 2 (w) 
is thus seen to have the same basic structure as the non- 
interacting g(ui). In particular, the u — > quasiparticle 
behaviour of the Green function is obtained by expanding 
S^w) about u = 0, 



G 2 (w) 



A R (w) 



(29) 



where e* = e; + |f7no + £ 2 (0) is the renormalized level 
energy, and 



Z 2 = 



1 - 



j=0J 



(30) 



is the quasiparticle weight. That G 2 (w) is a renormalized 
version of the non-interacting g(to) at low-frequencies is 
of course a direct reflection of adiabatic continuity to 
the non-interacting limit, in which sense the system is 
a generalized Fermi liquid (GFL). For some fixed ei and 
5, on switching on a small U, one expects the system to 
evolve smoothly into this perturbative state and as such 
to behave as a (local) Fermi liquid. As we now show, 
however, such behaviour does not arise at the particle- 
hole symmetric limit of the model. 



B. The particle- hole symmetric limit 

That the behaviour at the particle-hole symmetric 
point is qualitatively different from that described above, 
arises because (see sec.|TTJ the pole in the non-interacting 
d(u) lies precisely at loq = 0. On substituting eqn. (|15|) 
into eqn. ([TO)) , Im °II(a;) itself contains a pole at zero 
frequency: 



1 



-Im °n(w) = \Ql8(u) 



M 



+ Q d h (uj)+ du/d b (u/)d b (a/ - M). (31) 



Physically, this pole describes zero-energy spin-flip exci- 
tations of the impurity at the non-interacting level: the 



u> = pole in d(uj) reflects a state which is equally likely 
to contain a single ]- or J.-spin electron. 

The second-order self energy term £ 2 (w) is obtained as 
before by substituting eqn. (|3"Tj) into eqn. (|2"0|) . Whereas 
away from particle-hole symmetry the key result was a 
gap in S 2 (w) around the Fermi level (see eqn. (f2"3"]) ). here 
instead the pole in Im °U(ui) leads to a corresponding 
pole in the self-energy: 



S 2 M l "£ S kU 2 Ql 5(w). 



(32) 



Such behaviour is strikingly different from the GFL 
physics described above. The real part of T,^(uj) obtained 
from eqn. ((32)) diverges as oj — > 0, i.e. 



,o U 2 Ql 



(33) 



And the corresponding G 2 (w) = \g{oS) 1 — iUno — 
S 2 (w)] _1 is thus of form 



G 2 (w) 



Ioj I <<5 



1 



A R (w) 



(34) 



(noting that the renormalized level energy e* , given gen- 
erally by e* = 6i + E R (0), vanishes by symmetry at the 
particle- hole symmetric point t\ = —U/2). 

Eqn. I|34p clearly cannot be written as a renormalized 
version of the non-interacting g(u>), indicating that the 
particle-hole symmetric point is not perturbatively con- 
nected in U to the non-interacting limit. Instead, it is 
readily shown by expanding the hybridisation function 
Ar(oj) about oj — that 



G 2 (w) 



Qo 



(35) 



with a renormalized interaction U e s — UQq (and Qo = 
[1 + ^f-]^ 1 from eqn. HH)). The significance of this re- 
sult is that it is a renormalized version of the atomic limit 
(V = 0) propagator, Gal( w ) — [ z ~ \U 2 / z] -1 , indicative 
of the LM nature of the phase arising at particle-hole 
symmetry. Instead of the single pole inside the gap seen 
away from particle-hole symmetry, two poles thus arise, 
at ui = ± if/off, and representing as such renormalized 
versions of the atomic limit Hubbard satellites (them- 
selves occurring at uj = ±4f7). 

Even at the basic level of second-order perturbation 
theory, one sees then that introducing a gap in the con- 
duction band of the Anderson model changes significantly 
its low-energy physics, with non-analyticities arising in 
the perturbative self-energies in the particle-hole sym- 
metric limit. We now extend the treatment to arbitrary 
order in the interaction U, bolstering the conclusions 
drawn above and in doing so obtaining some exact re- 
sults for the problem. 



5 



t-O 



+ 




+ 




(a): E(fi))+A(ffl) 
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FIG. 3: Skeleton expansion of the self-energy. The exact 
is given by the sum of all skeleton self-energy diagrams 
constructed from the exact G(u). 



IV. INFINITE ORDER PERTURBATION 
THEORY IN U 

One cannot of course derive perturbatively the full fre- 
quency dependence of the impurity Green function to all 
orders, but its behaviour inside the gap can be ascer- 
tained exactly using some simple arguments. Our main 
result is to determine the conditions under which per- 
turbation theory in U is valid, and the resulting low- 
frequency behaviour of the impurity Green function. 

The argument follows that of the classic paper by 
Luttinger— , in which the skeleton expansion of the 
self-energy is used to obtain self-consistently the low- 
frequency behaviour of a wide class of interacting prob- 
lems to all orders in the interaction. The basis of the 
approach is that if perturbation theory holds, the exact 
self-energy is equal to the sum of all skeleton self-energy 
diagrams constructed from the exact Green function, as 
shown diagramatically in fig. [3] By analysing the general 
low-frequency behaviour of all such skeleton diagrams, 
we show that this condition is satisfied only if the exact 
Green function contains a single pole inside the gap at a 
non-zero frequency. 

A key insight of ref. HH is that the low- 
energy behaviour of the self-energy can be deduced 
perturbatively — via the skeleton expansion — from the 
low-energy single-particle excitations of the exact G(uj). 
To be more precise, we refer to the positive frequency ex- 
citations of D(w) as 'particle energies' and the negative 
frequency excitations as 'hole energies'. Then, it can be 
shown— quite generally that ImS(w) has weight at pos- 
itive frequencies corresponding to n + 1 particle energies 
minus n hole energies, or at negative frequencies corre- 
sponding to n + 1 hole energies minus n particle energies, 
with n > 1. Combining this result with the Dyson equa- 
tion [giving G(u>) in terms of E(w)] generates a pair of 
self-consistency equations for In what follows, we 

show that any self-energy that contains poles inside the 
gap (i.e. for \u>\ < 6) is inconsistent with these require- 
ments. 

Consider the situation in which contains at least 

one pole in the range —8 < lu < 8, and label by loq the 
frequency of the pole closest to uj = 0. In the following 
we take u)q > 0; similar reasoning can be applied if the 
closest pole is at a negative frequency. Then define u>\ 
to be the least negative frequency at which E J (w) has 
weight (either in the form of another pole, or the upper 
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(c): G(co) from (ii) 
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FIG. 4: Sketch showing how poles in E(cj) lead (via the Dyson 
equation) to poles in G(ui). (a) Solid lines represent the 
schematic form of E i (oj) + Ai(lj) with poles in the gap, dashed 
lines represent the corresponding E R (lj) + Ar(oj) from Hilbert 
transformation, and the dotted lines show the constructions 
used for determing the poles in G(ui) from Dyson's equation 
§5§. (b) The pole structure of G(ui) obtained from case (i) of 
(a), (c) The pole structure of G(u>) obtained from case (ii) of 
(a). 



edge of the negative-frequency continuum) , and u>2 to be 
the next frequency above luq at which S I (w) has weight. 
We sketch the resulting E : (w) + Ai(w) in fig-EJa), taking 
(without loss of generality) both lj% and u>2 to correspond 
to discrete poles rather than the edges of the continua. 

In fig. |U(a), we also sketch the E R (w) + Ar(w) corre- 
sponding to the £ : (o;) + At(o;) described above. Its form 
follows from the Hilbert transform eqn. ©. Within each 
frequency range in which S I (w) is zero, an argument akin 
to that before eqn. (|2rJ)) shows that E R (w) is monoton- 
ically decreasing. Moreover, at frequencies where S : (w) 
has poles, it is readily seen from eqn. © that E r (oj) 
diverges as 1/(lo — lo p ) (with lo p the particular pole fre- 
quency). Hence E R (w) + A(w) spans the full range from 
oo to -co between each excitation in E 1 ^) + Ai(w). 

The real and imaginary parts of the schematic self- 
energy can then be inserted into the Dyson equation to 
determine G(ui). Since Ai(w) > for \uj\ > 5, it follows 
from eqns. §5§ and §5§ that G((J) has poles only for 
\lo\ < 5, at frequencies given by the solution of u> — t\ = 

s R H + A R ,H. 

The resultant poles in G(uS) thus lie at the frequencies 
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where the dashed lines in fig. [4] intersect straight lines of 
unit slope, as shown by the dotted lines in fig.|U[a). From 
the diagram, it is clear that within the region uj\ < uj < 
u)2 there will always be two such poles in G(u>). Either 
there is one on each side of uj = 0, which we denote by 
'case (i)' and show schematically in fig. 0Jb); or both 
poles are at uj > 0, denoted 'case (ii)' (see fig. Etc)) [Note 
that if ujq = 0, then the only possible case is (i)]. 

Now consider the uj > behaviour of obtained 
by inserting this G{uj) into the skeleton expansion. As 
discussed above, this is obtained^ by considering the en- 
ergies of possible particle and hole excitations in G{uj). 
In case (i), the lowest particle energy is at an uj > ujq, 
and the lowest hole energy is at an u < 0. SinceiS the 
lowest-energy pole in for uj > is at an energy 

corresponding to two particle energies minus one hole en- 
ergy, this lowest-energy pole is necessarily at a frequency 
uj > Iloq and is therefore not consistent with the starting 
S I (o;) in fig. \^a). Likewise in case (ii), the lowest parti- 
cle energy is at an w > 0, and the lowest hole energy is at 
an uj < oj\. It follows that any resulting pole in S I (w) is 
at an uj > |ci>i| with |ci>i| > ujq by definition. Hence case 
(ii) is also not self-consistent. 

We have thus shown that there are no self-consistent 
perturbative solutions of the GAIM with a pole in S i (oj) 
at an ujq > 0. Similar arguments can be used to prove 
the same is true for a pole at ujq < 0. Therefore, any 
perturbative (Fermi liquid) solution cannot contain poles 
in £ ! (w) inside the gap. 

Finally we show that given the conclusion drawn 
above, particle-hole symmetry always leads to a break- 
down of perturbation theory. To that end, imagine start- 
ing with a perturbative S I (aj) [i.e. no poles in the gap] 
in the particle- hole symmetric limit. Its corresponding 
real part is monotonically decreasing inside the gap, and 
furthermore satisfies 6i + £ R (0) = by particle-hole sym- 
metry. Therefore the Green function obtained by substi- 
tuting this self-energy into the Dyson equation (|5|) must 
have a pole at w = 0. But then the S ! (w) obtained by 
substituting this G{uj) into the skeleton expansion would 
itself have a pole at zero frequency, breaking the pertur- 
bative self-consistency of the skeleton expansion. Hence 
the particle-hole symmetric point cannot be perturbatively 
connected to the non-interacting limit, and as such is a 
non Fermi-liquid, for all S > 0. 

Away from particle- hole symmetry by contrast, the 
pole in G(uj) arising from a perturbative S : (w) can (and 
will) be at a non-zero frequency. In this case the skele- 
ton expansion would not lead to a pole inside the gap in 
E r (o;), which is consistent with the form of the starting 
I] 1 (a;). Hence such perturbative solutions can in princi- 
ple exist; and the results of sec. HH naturally lead us to 
expect them, at least for some non-vanishing radius of 
convergence in U. 



V. THE LUTTINGER-WARD FUNCTIONAL 

By closing the skeleton self-energy diagrams in fig. 
[3] with an additional propagator line, one obtains the 
Luttinger-Ward functional $[C?] 21 , from which the exact 
self-energy follows from the functional derivative 



(36) 



g=g 



Using eqn. (|36|) . and that &[G] is invariant 21 under 
a frequency shift of all its propagators Q{uj), it can be 
shown that the 'Luttinger integral' Jl must vanish: 



h, = Im 



° d(J ^M G(w)= o. (37) 

-no OUJ 



This can be used in turn to derive a generalized Friedel 
sum rule 2 for the GAIM, which must be satisfied by 
any perturbative solution of the problem. Consider the 
change in the density of states upon introducing the im- 
purity, Ap(u>) = p(uj) - phost(^). This is given by 



Ap(uj) 



1 



Im 



^2G kk (uj)+G(uj)-J2 



1 



Z - £k 



(38) 



The equations of motion for the host Green functions 
generate 2 - 



k k 

= E 



i 



Z - 6k 



Z - 6k 



^— ' (Z - 6 k ) 2 



8A{uj) 
duj 



G(uj), 



and the Dyson equation ([5]) gives 



d_ 



Giuj)- 1 = 1 



(39) 
(40) 

(41) 



Combining eqns. 
result^ 



duj duj 
(|4"Tj) thus leads to the general 



Ap{uj) 



1 



Im 



dlnGH" 1 dE(u) 



duj 



du> 



G(uj) 



(42) 



The number of electrons introduced by the impurity, 
iii mp , is now obtained by integrating 2Ap(uj) up to uj = 
(the factor of 2 arising from a spin sum). From eqn. (|41[) . 
and using the Luttinger integral condition eqn. (f37|) , one 
has 



Im 



8 

duj— lnG(w)- 1 . (43) 



From equation ((5]) it follows that 

'uj - a - A r (uj) - S R (w) 



ImlnG(w) 



7T 

tan 

2 



Ai(w) + S : (w) +n 



(44) 



7 



(taking the principal range — 7r/2 < tan 1 x < 7r/2), and 
hence^ 



= 1 tan 



£i + A R (0) + S R (0) 
Ai(0) + S I (0) + ?7 



(45) 



This Friedel sum rule is quite general for the Ander- 
son impurity model. For the GAIM in particular one has 
Ai(0) = = Ar(0), and the results for the self-energy 
deduced in the previous section showed that E 1 (0) neces- 
sarily vanishes if perturbation theory holds (GFL phase). 



Hence n- u 



— tan 1 



-E R (0) 



from which 



ip = 1 - sgn(e*) 



(46) 



with e* = 6j + E R (0) the renormalized level energy. For 
any e* < 0, the additional number of electrons induced 



by the impurity is thus 



2. while riy 



for any 



e* > 0. This behaviour — exclusively integral values of 
-is physically natural for a system with a gapped 



spectrum, although it is of course in marked contrast to 
the metallic (gapless) Anderson model for which Ai (0) > 
and hence (from eqn. (|45p ) ni mp is a continuous function 
ofe*. 



VI. CONCLUSION 

Finally, we tie together the results of the previous sec- 
tions. We have shown that self-consistent perturbative 
solutions of the problem must necessarily have a non-zero 
renormalized level, e*, since this leads to a single pole in 



G{u>) at a non-zero frequency and hence a perturbatively- 
consistent S(w) from the skeleton expansion. If e* is 
negative the pole in G(ui) is at a negative frequency and 
"imp = 2; while if e* is positive the pole is at a positive 
frequency and n- lmp = 0. 

In sec. HVi we deduced that such perturbative solutions 
of the problem can only exist away from particle-hole 
symmetry: it is impossible to construct a G(u>) that is 
both consistent with particle-hole symmetry and the per- 
turbative skeleton expansion of the self-energy. Hence the 
particle-hole-symmetric point of the problem is a non- 
Fermi-liquid LM phase, and any approximation to the 
GAIM must therefore take this into account 1 -. Away 
from particle-hole symmetry, there is nothing to prevent 
the existence of perturbative GFL solutions, and these 
are indeed found numerically in e.g. the NRG^. 

Finally we note that since the derivation of the 
generalized Friedel sum rule, eqn. (|46|) above, has 
been obtained from a perturbative construction of the 
Luttinger Ward functional it is not valid in the 

LM phase. At the particle-hole-symmetric point, it is 
readily shown instead that n; mp = 1 by symmetry. In 
our forthcoming local moment approach to the models, 
we argue that n lmp = 1 is a general result that also holds 
for the LM phase away from particle-hole symmetry, 
and present results of this non-perturbative approach 
for the dynamics of the problem in the two phases. 

Acknowledgments 

We are grateful to the EPSRC for financial support. 



1 P. W. Anderson, Phys. Rev. 124, 41 (1961). 

2 A. C. Hewson, The Kondo Problem to Heavy Fermions 
(Cambridge University Press, Cambridge, 1993). 

3 T. K. Ng and P. A. Lee, Phys. Rev. Lett. 61, 1768 (1988). 

4 L. I. Glazman and M. E. Raikh, Sov. Phys. JETP Lett. 
47, 452 (1988). 

5 D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, 
D. Abusch-Magder, U. Meirav, and M. A. Kastner, Na- 
ture 391, 156 (1998). 

6 S. M. Cronenwett, T. H. Oosterkamp, and L. P. Kouwen- 
hoven, Science 281, 540 (1998). 

7 L. P. Kouwenhoven et al., in Mesoscopic Electron Trans- 
port, edited by L. L. Sohn (Kluwer, Dordrecht, 1997). 

8 D. L. Cox and A. Zawadowski, Advances in Physics 47, 
599 (1998). 

9 R. Bulla, T. Costi, and T. Pruschke, cond-mat/0701105 
(2007). 

10 O. Sakai, Y. Shimitzu, H. Shiba, and K. Satori, J. Phys. 
Soc. Jpn. 62, 3181 (1993). 

11 A. V. Balatsky, I. Vekhter, and J.-X. Zhu, Rev. Mod. Phys. 



78, 373 (2006). 

12 J. Bauer, A. Oguri, and A. C. Hewson, J. Phys.: Condens. 
Matter 19, 486211 (2007). 

13 T. Saso, J. Phys. Soc. Jpn. 61, 3439 (1992). 

14 J. Ogura and T. Saso, J. Phys. Soc. Jpn. 62, 4364 (1993). 

15 C. C. Yu and M. Guerrero, Phys. Rev. B 54, 8556 (1996). 

16 K. Takegahara, Y. Shimitzu, and O. Sakai, J. Phys. Soc. 
Jpn. 61, 3443 (1992). 

17 K. Chen and C. Jayaprakash, Phys. Rev. B 57, 5225 
(1998). 

18 J. M. Luttinger, Phys. Rev. 121, 942 (1961). 

19 M. R. Galpin and D. E. Logan, Eur. Phys. J. B 62, 129 
(2008). 

20 M. T. Glossop and D. E. Logan, J. Phys.: Condens. Matter 
14, 6737 (2002). 

21 A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinski, 
Methods of Quantum Field Theory in Statistical Physics 
(Dover, New York, 1977). 



